Scaling of Entanglement Entropy in the Random Singlet Phase 
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We present numerical evidences for the logarithmic scaling of the entanglement entropy in critical 
random spin chains. Very large scale exact diagonalizations performed at the critical XX point 
up to L = 2000 spins i lead to a perfect agreement with recent real-space renormalization-group 
predictions of Refael and Moore [Phys. Rev. Lett. 93, 260602 (2004)] for the logarithmic scaling of 
the entanglement entropy in the Random Singlet Phase with an effective central charge c = c x In 2. 
Moreover we provide the first visual proof of the existence the Random Singlet Phase with the help 
of quantum entanglement concept. 



PACS numbers: 75.10.Pq, 03.67.Mn, 75.10.Nr 



The study of quantum phase transitions through quan- 
tum entanglement concepts provides a new way to un- 
derstand strongly correlated systems near criticality. In 
one dimensional systems, such as quantum spin chains, 
entanglement estimators exhibit universal features close 
to a critical pointi*^. One of this estimator is the en- 
tanglement entropy of a subsystem A with respect to a 
subsystem B. Defined as the Von Neumann entropy of 
the reduced density matrix for either subsystem 



S = —TrpA In Pa = -Trp B \np B , 



(1) 



this quantity displays very interesting scaling behavior 
for conformally invariant critical theories in one dimen- 
sion (ID). Indeed, as shown first by Holzey, Larsen and 
Wilczek^ in the context of geometric entropy related to 
black hole physics, the entanglement entropy of a sub- 
system of length x embedded in an infinite system is ex- 
pected to scale like 



S(x) ~ — lax. 



(2) 



The number c is the so-called central charge which is, 
for instance, for the critical XXZ spin-^ chain cxxz = 1 
or for the spin-i Ising chain in transverse field at criti- 
cality ci s i n g = 1/2. This result [Eq. J5J] has been veri- 
fied numerically^ as well as analytically in Ref£ where 
some simple connections have been established between 
thermodynamic entropy and entanglement entropy. An 
important extension to critical and non-critical systems 
with finite size, finite temperature and different boundary 
conditions has been achieved by Calabrese and Cardy^. 
They showed for instance that for critical systems of fi- 
nite size L with periodic boundary conditions, Eq. i|2Jl 
should be replaced by 



c / L . ttx , 
o{L,x) = - In -sm(— J 

O \ 7T L 



•si, 



(3) 



where si is a constant related to the UV cut-off. 

Although such a logarithmic scaling of the entangle- 
ment entropy seems closely related to the conformal in- 
variance of the critical system, it has been shown recently 
by Refael and Moored that such a critical scaling is also 



expected for some random critical points. Indeed, using 
an analytic real-space renormalization-group (RSRG) ap- 
proach, they have shown that random critical spin chains 
display similar features that clean ones with an effective 
central charge c — c x In 2 so that 



S(x) 



In a; 4- constant. 



(4) 



This surprising results can be derived using the RSRG 
method introduced by Ma, Dasgupta and Hu£ several 
years ago to study random spin chains. As a result, 
any amount of randomness introduced as a perturba- 
tion in a clean XXZ critical spin-i chain is relevan-feiS 
and drives the system to the so-called Random Singlet 
Phase (RSP)£, associated with an infinite randomness 
fixed point (IRFP) for the RSRG transformation 9 . The 
RSP can be depicted as a collection of singlet bonds of 
arbitrary length (see Fig. Then, utilizing the very 





FIG. 1: Schematic picture for the entanglement entropy of a 
subsystem of length x in the random singlet phase. The entan- 
glement is just due to the singlets connecting the subsystem 
with the rest of the chain. In this picture, S(x) — 5 x In 2. 



simple result that the entanglement entropy of a spin ^ 
involved in a singlet with its partner is In 2, in the RSP 
the entanglement of a segment with the rest of the system 
is just given by In 2 times the number of singlets which 
cross the boundary of the segment, as depicted in Fig.^ 
Using this fact as well as an accurate RSRG calculation, 
Refael and Moore have then been able to determine pre- 
cisely that the number of singlets connecting a segment 
of size x with the rest of the system is (1/3) lnx, leading 
to the formula Q . The purpose of this communication is 
to investigate numerically the entanglement of the RSP 
and compare the RSRG prediction [Eq. |@J] with exact 
computations. 
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Exact computation of the entanglement entropy — In 
order to compute the entanglement entropy of a subsys- 
tem, one needs to calculate the corresponding reduced 
density matrix. For general XXZ spin chains governed 
by 



WxXZ - J ^2 O ( S -t S 0+l + S 3 S f+l) 



ASjSj +1 



(5) 

the non-critical regime (achieved if |A| > 1) can be in- 
vestigated using the corner transfer matrices of the cor- 
responding two-dimensional (2D) classical problemii*i£. 
On the other hand, along the critical line (—1 < A < 1), 
an analytical computation of S(x) is more difficult and 
conformal field theory (CFT) tools are then required 6 . 
Another alternative consists in performing numerical ex- 
act diagonalizations (ED) of finite lengths spin chains, 
but it is limited to L max ~ 40 spins i when A ^ 
Nevertheless, the XX point A = is special because 
the spin Hamiltonian can be rewritten using the Jordan- 
Wigner transformation as a free-fermions model 



Ti-xx = 



J 



£ 

3 



[4 



(6) 



for which the density matrix can be expressed as the 
exponential of a free-fermion operator*^. It turns out 
that the reduced density matrix is completely determined 
by the x x x correlation matrix C(x) 7 defined by 



((c\c x ) 



C(x) = 



(c\c 2 ) 



(c|ci) 



(4 



2 C V 



V 



{c\c x )\ 



(7) 



The matrix elements Cy = ( c l c j) can be calculated either 
numerically by diagonalizing the free-fermion Hamilto- 
nian in momentum space or analytically in some special 
casesi^. The entanglement entropy of a subsystem of size 
x embedded in a larger system is then given by 



S(x) = -J2 [A* In A fe + (1 - A fe ) ln(l - A*)] . 



(8) 



where the Xk are the eigenvalues of C(x). 

Let us now concentrate on the disordered XX spin-^ 
chain, governed by the random hopping Hamiltonian on 
a periodic ring of length L 



T~lxx 



£4 



1 + C j+l C 3 



[c\ Cj 

3=1 

+Jl exp(z7rAT)(c| j ci + c\cl) 



(9) 



where Jj are positive random numbers chosen in a flat 
uniform distribution within the interval [0, 1}^, and the 
second term in the right hand side ensures that peri- 
odic boundary conditions are imposed in the spin prob- 
lem. The total number of fermions is M — L/ 2 in 



the ground-state (GS). The way to diagonalize Hxx is 
straightforward and has already been explained by sev- 
eral authorsiSiiS. As a check, we have first computed the 
entanglement entropy (JSJ for clean systems (i.e. Ji is a 
constant) of total sizes L = 500 and L = 2000. Tech- 
nically, this only involves computing the elements (c\cj) 
by diagonalizing the free-fermions Hamiltonian ©, and 
then one needs to diagonalize C [Eq. J7J] using standard 
linear algebra routines^. The results are shown in Fig. [21 
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FIG. 2: (color online) Entanglement entropy of a subsystem 
of size x embedded in a closed ring of size L, shown vs x in a 
log-linear plot. Numerical results obtained by exact diagonal- 
izations performed at the XX point. For clean non random 
systems with L = 500 and L — 2000 (open circles), S(x) is 
perfectly described by Eq. Q (red and blue curves). The data 
for random systems have been averaged over 10 4 samples for 
L = 500, 1000, 2000 and 2 x 10 4 samples for 100 < L < 400. 
The expression 0.8595 + =|2 In a; (dashed line) fits the data in 
the regime where finite size effects are absent. 



where we can see that S(L,x) is perfectly described by 
the CFT prediction Eq. 0. Note also that the constant 
term is found to be si ~ 0.726, in excellent agreement 
with the recent analytical prediction of Jin and Kore- 
pm*i. 

For the random case, the same technique has been used 
but a bigger computational effort was necessary to aver- 
age over a large number of independent random samples. 
Practically the number of samples used was 2 x 10 4 for 
L = 100, 200, 300, 400 and 10 4 for L = 500, 1000, 2000 
which required 2000 hours of CPU computational time. 
The results for the disorder averaged entanglement en- 
tropy are shown in Fig. |2] When the subsystem size x 
is large enough (typically x > 20), the expression (gj 
derived by Refael and Moore describes perfectly the be- 
havior of the disorder average entanglement entropy, i.e. 
a logarithmic scaling with an effective central charge 
c = In 2. One can notice that when the subsystem size 
approaches L/2 some finite size effects are visible, as it 
is the case in clean systems. 
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Signature of the Random Singlet Phase — The very 
good agreement found between exact numerical diago- 
nalizations and RSRG calculations for the entanglement 
properties in the RSP is a new proof in favor of the ran- 
dom singlet nature of the GS, also supported by recent 
neutron scattering experiments performed on the dis- 
ordered spin chain compound BaC^^io.sGeo.s^O-^. 
Another way to get more insight on these long distance 
effective singlets in the GS consists in looking at the prob- 
ability distribution of the entanglement entropy. Indeed, 
since each singlet is expected to contribute as a In 2 in 
the entanglement entropy, we can focus on the probabil- 
ity distribution of S / In 2 for a given subsystem embedded 
in a larger system. In order to get a correct statistical pic- 
ture for the typical behavior of this random singlets for- 
mation, one needs a huge number of disordered samples. 
We chose to study 10 5 independent realizations. The 
price to pay is that not too large systems can then be di- 
agonalized. Nevertheless, only focusing on L = 100 spins 
is enough to get good insights on the RSP. Indeed, in- 
stead of increasing the system size to achieve the physics 
of the RSP, according to the disorder induced crossover 
phenomena observed for the RSRG flowii one can rather 
keep L fixed and increase the disorder strength to get 
closer to the IRFP and therefore deeper in the RSP. Let 
us thus consider strong disorder distributions for the cou- 
plings Ji, like 

v(j) = - 5 j- 1+s '\ (10) 

parametrized by a disorder strength 5 > 1. This distri- 
bution is quite natural to mimic strong disorder effects 
since at the IRFP, the fixed point distribution for the 
random couplings is achieved for <5 — > oo. 

I order to minimize the finite size effects, we consider 
half of the chain as a subsystem and compute S(L, L/2) 
for each sample. Nevertheless, in order to get a good 
understanding, it is important to notice that the par- 
ity of L/2 is crucial. Indeed if L/2 is odd only an odd 
number of random singlets can connect both subsystems 
whereas if L/2 is even, the number of cut singlets will be 
even, none singlet being also a possibility. This fact is 
actually clearly visible in Figs. where we have plotted 
the probability distributions P evon = P(S(100, 50)/ In 2) 
(green histogram) as well as P odd = P(S(100, 49)/ In 2) 
(red histogram), for 6 = 1, 2.5, 5, 10. Whereas for 
6 = 1 (Fig. Eta)) P ovon (P odd ) displays an integer-peaks 
structure, signature of the RSP, only for S/ In 2 = 2 
(<S/ln2 = 3) and that a non-negligible statistical weight 
lies between for non-integer values, when the disorder 
strength increases, the integer-peaks structure becomes 
more and more pronounced as visible in Figs. |2Ib-d). 
The combined distributions p cvcn -f p°dd are a i so plot- 
ted in the insets of Figs. El Thanks to the entanglement 
entropy, we provide a clear visual proof for the RSP. 

Discussion and conclusion — Non disordered critical 
spin chains can be described by a conformally invariant 
field theory from which an universal number c, the central 
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FIG. 3: (color online) Probability distribution for the entan- 
glement entropy in the random singlet phase for the disor- 
dered XX spin-i chain obtained by exact diagonalizations on 
L — 100 spin chains. Green histograms correspond to an even 
subsystem size with 50 sites (P cvcn = P(<S(100, 50)/ In 2)) and 
red histograms correspond to an odd subsystem with 49 sites 
(podd = p(s(ioo,49)/ln2)). For each' disorder strength, 
S = 1 (a), 8 = 2.5 (b), 5 = 5 (c), 5 = 10 (d), 10 5 differ- 
ent random samples have been diagonalized. The insets show 
the combined distributions P cvcn + P odd . Note that all the 
distribution functions have been normalized to unity. 



charge, emerges. This central charge, also called confor- 
mal anomaly number, appears in the leading finite size 
(or finite temperature) correction to the free energy^ as 
well as in the entanglement entropy ©. The power-law 
behavior for the spin-spin correlations functions is also 
universal with well defined critical exponents^ as well as 
exact amplitude^. 

In the case of random critical chains, while the RSRG 
framework provides universal critical exponents, the am- 
plitudes of correlations functions are non-universal num- 
bersSi. On the other hand, the RSRG treatment for 
the entanglement entropy provides the exact prefactor 
equal to This prediction has been checked numer- 
ically using exact numerical diagonalizations on large 
scale random critical spin chains. The perfect agreement 
between exact simulations and the perturbative RSRG 
provides, to the best of our knowledge, the first example 
of an exact critical amplitude computed within this tech- 
nique. It is also interesting to notice that this finding 
of c = In 2 < 1 would be consistent with a generalized 
c— theorem built on entanglement concepts for non con- 
formal random critical points^. Nevertheless, the identi- 
ficaton of other physical quantities besides entanglement 
that are controlled by this number c turns out to be very 
challenging and more subtle than using a simple anal- 
ogy with the clean case. Indeed, since in conformally 
invariant clean systems a non universal velocity factor 
appears tied to c in the usual thermodynamic quantities 
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like the specific heafj^ 3 . or the aforementioned correction 
to the free energy, the analogy here breaks down because 
the velocity of excitations is not defined anymore in the 
RSP. 

To conclude, we believe that the results presented 
in this communication provide a new insight on the 
random singlet phase as well as a first visual proof of 
the large scale effective singlets formation. The fact 
that even at the infinite randomness fixed point the 
entanglement entropy still scales logarithmically with 



the subsystem size provides a non trivial extension of the 
quantum entanglement concepts to random quantum 
critical points. 
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